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Abstract. As an alternative explanation of the dimming of distant supernovae it has 
recently been advocated that we live in a special place in the Universe near the centre 
of a large void described by a Lemaitre-Tolman-Bondi (LTB) metric. The Universe 
is no longer homogeneous and isotropic and the apparent late time acceleration is 
actually a consequence of spatial gradients in the metric. If we did not live close 
^ ' to the centre of the void, we would have observed a Cosmic Microwave Background 

C^ . (CMB) dipole much larger than that allowed by observations. Hence, until now it has 

been argued, for the model to be consistent with observations, that by coincidence 

we happen to live very close to the centre of the void or we are moving towards it. 

'^ , However, even if we are at the centre of the void, we can observe distant galaxy clusters, 

pg ' which are off-centre. In their frame of reference there should be a large CMB dipole, 

CO ' which manifests itself observationally for us as a kinematic Sunyaev-Zeldovich (kSZ) 

effect. kSZ observations give far stronger constraints on the LTB model compared to 
C^^ ' other observational probes such as Type la Supernovae, the CMB, and baryon acoustic 

^^ , oscillations. We show that current observations of only 9 clusters with large error bars 

f^ I already rule out LTB models with void sizes greater than ~ 1.5 Gpc and a significant 

underdensity, and that near future kSZ surveys like the Atacama Cosmology Telescope 
(ACT), South Pole Telescope (SPT), APEX telescope, or the Planck satellite will be 
/\ ' able to strongly rule out or confirm LTB models with giga parsec sized voids. On the 

jrt ' other hand, if the LTB model is confirmed by observations, a kSZ survey gives a unique 

possibility of directly reconstructing the expansion rate and underdensity profile of the 
void. 
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1. Introduction 

Distant supernovae appear dimmer than expected in a matter-dominated homogeneous 
and isotropic FRW universe. The currently favoured explanation of this dimming is 
the late time acceleration of the universe due to a mysterious energy component that 
acts like a repulsive force. The nature of the so-called Dark Energy responsible for the 
apparent acceleration is completely unknown. Observations seem to suggest that it is 
similar to Einstein's cosmological constant, but there is inconclusive evidence. There 
has been a tremendous effort in the last few years to try to pin down deviations from a 
cosmological constant, e.g. with deep galaxy catalogues like 2dFGRS P and SDSS [2], 
and extensive supernovae surveys like ESSENCE [3], SNLS [1], and SDSS-SN [5], and 
many more are planned for the near future e.g. DES [B], PAU [3, E], BOSS [S] and 

JDEM [in]. 

In the meantime, our realisation that the universe around us is far from 
homogeneous and isotropic has triggered the study of alternatives to this mysterious 
energy. Since the end of the nineties it has been suggested by various groups [HI [T21 
[T3l [HI [15l dni [TTl dH dm that an isotropic but inhomogeneous Lemaitre-Tolman-Bondi 
universe could also induce an apparent dimming of the light of distant supernovae, in 
this case due to local spatial gradients in the expansion rate and matter density, rather 
than due to late acceleration. It is just a matter of interpretation which mechanism is 
responsible for the dimming of the light we receive from those supernovae. Certainly the 
homogeneous and isotropic FRW model is more appealing from a philosophical point 
of view, but so was the static universe and we had to abandon it when the recession of 
galaxies was discovered at the beginning of last century. 

There is nothing wrong or inconsistent with the possibility that we live close to 
the centre of a giga parsec scale void. Such a void may indeed have been observed as 
the CMB cold spot [201 [211 [22] and smaller voids have been seen in the local galaxy 
distribution [23,121]. The size and depth of the distant observed voids, i.e. tq ^ 2 Gpc 
and Qm ~ 0.2 within a fiat Einstein-de Sitter universe, seems to be consistent with that 
in which we may happen to live [25], and could account for the supernovae dimming, 
together with the observed baryon acoustic oscillations and CMB acoustic peaks, the 
age of the universe, local rate of expansion, etc. [T9] . 

Moreover, according to the theory of eternal inflation [22], rare fluctuations at 
the Planck boundary may be responsible for the non-perturbative amplification of 
local inhomogeneities in the metric, which would look like local voids in the matter 
distribution [27]. In the eternal inflation approach one assumes to be a typical observer 
whose local patch comes directly from a rare fluctuation within an inflationary domain at 
the Planck scale. Moreover, since it is a rare fluctuation it should be highly spherically 
symmetric, and the theory predicts that we should live close to the centre of such a 
void [27]. The size and depth of those voids depends on the theory of inflation on very 
large scales and may be a probe (perhaps the only probe) of the global structure of the 
universe. 
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In fact, observations suggest that if there is such a large void, we should live close 
to the centre, otherwise our anisotropic position in the void would be seen as a large 
dipole in the CMB. Of course, we do observe a dipole, but it is normally assumed to 
be due to the combined gravitational pull of the Virgo cluster, and the Shapley super 
cluster. There is always the possibility that we live off-centre and we are moving towards 
the centre of the void, so that the two effects are partially cancelled, giving rise to the 
observed dipole. However, such a coincidence could not happen for all galaxies in the 
void and, in general, clusters that are off-centred should see, in their frame of reference, 
a large CMB dipole. Such a dipole would manifest itself observationally for us as an 
apparent kinematic Sunyaev-Zeldovich effect for the given cluster. 

It is the purpose of this paper to study the very strong constraints that present 
observations of the kSZ effect already put on the LTB void models, and predict how 
near future observations from kSZ surveys like ACT or SPT will strongly rule out (or 
confirm) LTB models with giga parsec sized voids. 

In section 2 we describe the general LTB void models, giving the corresponding 
Einstein-Friedmann equations, as well as parameterisations of their solutions. In a 
subsection we describe the GBH constrained model, where we assume the Big Bang 
is homogenous, and thus the model depends on a single function, the inhomogeneous 
matter ratio VtMij)- In section 3 we study the induced dipole for off-centred clusters and 
compute the size of the analogue velocity of those clusters depending on the parameters 
of the void model. In section 4 we analyse present observations and give constraints on 
the model from current observations. In section 5 we then explore the prospects that 
future experiments like ACT, SPT and Planck will provide for strongly constraining or 
even ruling out LTB models of the universe. Finally, in section 6 we give our conclusions. 

2. Lemaitre-Tolman-Bondi void models 

The Lemaitre-Tolman-Bondi model describes general radially symmetric space-times 
and can be used as a toy model for describing voids in the universe. The metric is 

ds^ = -dt^ + X^ (r, t) dr'^ + A^ (r, t)dn^ , (1) 

where dVl? = dO"^ + siv? 9d(j)'^ . Assuming a spherically symmetric matter source with 
negligible pressure, 

T!: = -pM{r,t)Si;5l (2) 

the (0,r) component of the Einstein equations, G*° = 0, implies X{r,t) = 
A'(r,t)/ ^yl — k(r), with an arbitrary function k(r) playing the role of the spatial 
curvature parameter. The other components of the Einstein equations read [281 [T6l [T9] 

A^ + k AA' k'ir) 

A^ + 2AA + k{r) = . (4) 
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Integrating the last equation, we get 

^2 ^^^) ^^^ 



(5) 



^2 ^3 j^2 

with another arbitrary function F{r), playing the role of effective matter content, which 
substituted into the first equation gives 
F'(r) 



Svr GpM(r,t) . 



A'A^{r,t) 
We can also use Eq. ([3]) to define the critical density as 






AA' 

'AA' 



ItxG pc{r,t). 



(6) 



(7) 



The boundary condition functions F{r) and k{r) are specified by the nature of the 
inhomogeneities through the local Hubble rate, the local total energy density and the 
local spatial curvature, 

A{r,t) 



H{r,t) 
F{r) = 
k{r) 



A(r,t)' 
Hl{r)nM{r)Al{r) 



H^{r)(nM{r) 



1 ) Alir) 



(8) 

(9) 

(10) 



where functions with subscripts correspond to present day values, Aq^t) = A{r,tQ) 
and Hoi^r) = H{r,to). With these definitions, the r-dependent Hubble rate is written 
as EHlIlB] 



H\r,t) = H',{r) 



^Aiir) 



Aa(r) 



^Mir)) 



Ao{r) 



A{r,t)J ' '"' " \A{r,t) 

and we fix the gauge by setting A^ir) = r. 

For light travelling along radial null geodesies, ds"^ = dVl? = 0, we have 
dt _ A'{r,t) 
dr 



(11) 



T- 



(12) 



(13) 



Vl - k{r) 
which, together with the redshift equation, 

c;iog(l + ^) _ A'{r,t) 
dr ~ ^1 - k{r) 

can be written as a parametric set of differential equations, with TV = log(l + z) being 
the effective number of e-folds before the present time, 
dt _ A'{r,t) 

dN~~A'{r,t) ' 

dr -y/l — k{r) 

dN~ A'{r,t) 



(14) 
(15) 
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2.1. The constrained GBH model 

In general LTB models are uniquely specified by the two functions k{r) and F{r) or 
equivalently by Ho{r) and r2j\/(r), but to test them against data we have to parameterise 
the functions, to reduce the degrees of freedom to a discrete set of parameters. For 
simplicity in this paper we will use the constrained GBH model [19] to describe the void 
profile. First of all, it uses a minimum set of parameters to make a realistic void profile, 
and secondly, it is assumed that the time to the Big Bang is constant for spatial slices. 
The second condition gives a relation between Ho{r) and QM{f), and hence constrain 
the models to one free function, and a proportionality constant describing the overall 
expansion rate. Our chosen model is thus given by 

n,dr) = n... + ( 1^. - i^o. ) ( ' '-'^''!'Kr'y'^'^ ) (16) 



m(r) = Ha 




2[fi/<(r)]^ 



(17) 



_fiA'(r) ,/9^ ^nM{r)\ "^(2n+l)(2n + 3) 

where flxif) = 1 ^ ^m{i^), and the second equation follows from the requirement of a 
constant time to a homogeneous Big Bang. We use an "inflationary prior" , and assume 
that space is asymptotically fiat, i.e. in the following we set fiout = 1- The model has 
then only four free parameters: The overall expansion rate Hq, the underdensity at the 
centre of the void Q-^, the size of the void ro, and the transition width of the void profile 
Ar. For more details on the model see Ref. 



3. The CMB sky seen by off-centre observers 

Imagine a cluster well embedded in a big void, but not exactly at the centre, and 
consider photons reaching the cluster from the last scattering surface (LSS). Because of 
the symmetry of the problem the smallest and largest redshifts are found along the radial 
direction. There are two effects contributing to the redshift of photons passing through 
the void: The dominant effect is caused by the higher expansion rate inside the void (see 
Eq. (IT4l) ). The photons coming from the farthest end of the void, crossing the centre, are 
inside the void for the longest time, and thus have the biggest redshift, while the photons 
arriving directly from the LSS along a radial geodesic will be affected for the shortest 
amount of time, and consequently suffer the least redshift. There is also a subdominant 
effect due to the change in the gravitational potential or the matter density: Photons 
coming from the farthest end of the void are first gravitationally redshifted when entering 
the void, and then gravitationally blueshifted after crossing the centre. On the other 
hand, photons arriving directly from infinity are only redshifted. There is a difference 
in the two redshifts, because of the time dependence of the underdensity, and hence 
the gravitational potential, with the subdominating effect leading to a small blueshift 
towards the centre. This effect is in some sense a large scale Rees-Sciama effect [29], see 

Fig.m 
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Figure 1. An off-ccntrc cluster of galaxies in a void will "observe" CMB photons 
coming from the last scattering surface from all directions. Due to the higher expansion 
rate inside the void, photons arriving through the centre (from the right in the figure) 
will have a larger redshift (Azin), than photons arriving directly from the LSS (left, 
with Azout)- There is a subdominant effect due to the time-dependent density profile 
(the solid line corresponds to the current time, while the dot-dashed line to one tenth of 
the present time). With a larger underdensity at later times, we have Azi > Az4, and 
Az2 -t- Az3 < 0, giving an overall difference Azi > Az2 + Az3 + Az4 or, equivalently, a 
subdominant dipole with a blueshift towards the centre of the void. The overall effect 
is a blueshift away from the centre. 



Consequently, in the ideal case of a spherical void, and a well embedded cluster, the 
cluster observer will see an almost perfect dipole in the CMB, aligned along the radial 
direction, and with the blueshift pointing away from the centre of the void, where the 
observer is (see Fig. [T]). The detailed effect of a spherical void on the CMB sky of an 
off-centre observer has been calculated in [30], and it is shown that for small distances 
from the centre of the void the dominating term is a dipole. For simplicity, in this paper 
we will estimate the change on the CMB sky as a pure dipole. To find the amplitude it 
is then enough to integrate two radial light rays one going towards the centre, and the 
other away from the centre of the void. The calculation of the grid of parameters is done 
using our public easyLTB program [31], but we have also used the completely different 
approach of Taylor expansion around an Einstein de Sitter solution, as detailed in |19j . 
and checked that we get the same result for the parameters in Fig. [21 It should be noted 
that this dipole approximation breaks down when the effective size of the void on the 
sky, as observed from the cluster, becomes too small (i.e. less than ~ 27r), or possibly 
when the density or Hubble expansion profiles have very contrived time dependence. 

Motivated by the observed size of our CMB dipole, we assume that we are located 
at the centre of the void, and observe clusters in the light cone at different redshifts. 
Each cluster is then observed at a certain time tci{z), related to its redshift. In order 
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to find the dipole as seen by the cluster, we integrate Eq. (IT4l) along the radial axis 
in the positive and negative direction, and always backwards in time, with a starting 
time tci{z), and an ending time tiss ~ 10^ yr. The size of the dipole can be easily 
calculated from the temperature seen by observers in the cluster in different directions, 
T{6) = T*/(l + z(0)), where T* is the temperature of the LSS. If we now look in opposite 
directions, towards and away from the centre of the void, we find 

AT _ T{9) -T _ [Tin -Tout I _ |^in - ^outl ,, gs 

T dipole T -Mn + ^out 2 + ^jn + ^out 

where T = (Tin+Tout)/2 is the mean temperature observed at the location of the cluster, 
and 2in/out are the redshifts to the LSS for radially ingoing/outgoing light rays. 

3.1. The kinematic Sunyaev-Zeldovich effect 

The hot gas (mainly the electrons) inside a cluster inverse Compton scatters CMB 
photons, changing their frequency distribution. The thermal Sunyaev-Zeldovich 
effect [32l[33] is the main effect, redistributing low energy photons to higher energies due 
to upscattering by the hot cluster gas, but also the intra-cluster gas works as a mirror 
rescattering photons from all directions towards the observer. If the CMB sky observed 
by the cluster is different than the CMB sky observed by us, there will be an additional 
change in the spectrum, which is known as the kinematic Sunyaev-Zeldovich effect |34j . 
Because of rotational symmetry around the axis associated with the line of sight, only 
changes projected along this axis will have an observational impact. In the LTB model, 
because the observer is supposed to be at the centre, an off-centred cluster will observe 
a CMB dipole exactly along the line of sight, which to first order in the perturbation 
is indistinguishable from a kinematic dipole due to the peculiar velocity of the cluster. 
The change in the photon intensity is given as 



and thus 



^ = -^ farnj^-dl, (19) 

4 e^ - 1 J c 



Trad (e^ - 1)^ ./ c (e^ - 1 



- CTTUe— ■ dl ^ -fjp Te J-;. z^ , (20) 



where x = hu/kT, t^ is the cluster optical depth, and the apparent peculiar velocity is 
related to the temperature dipole as 

/3, = ^ = ^ . (21) 

C 1 dipole 

Even though the kSZ signal from a single cluster can be interpreted as a peculiar velocity, 
a large void would result in a systematic trend for all clusters (with a certain scatter given 
by the intrinsic peculiar velocities of the clusters), with an average redshift or positive 
apparent velocity (see Fig. [2]). The average velocity as a function of the distance would 
give a direct and unique handle on the density and velocity profile of the void. At the 
same time, the absence of such a systematic average peculiar velocity in future kSZ 
surveys can strongly rule out a void of any appreciable size. 
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n,„=0.1 3, ro=2.D, Ho=0.65, Ar/ro=0.1 -D.9 (bottom to top in rigtit side) ii,„=0.13-0,32 (top to bottom), ro=2, Ho=0.65, Ar/ro=0.5 



ii,^=0.13, ro=1 .5-2.5 (top to bottom). Ho=0.65. Ar/ro=0.5 




Figure 2. Examples of the size of the dipole for different parameters of the constrained 
GBH model [11]. The dashed line in the left figure is the first order approximation 
given in [50] , 




fij„=0.23. ro=1.8, H„=0.65, Ar/r„=0.35 




Figure 3. The angular and redshift distribution of current observations together with 
the predicted dipole distribution for a giga parsec sized void model. Red triangles and 
blue squares represent positive and negative peculiar velocities respectively, with the 
size of the symbol indicating the magnitude of the velocities. 



4. Constraints from current observations 



We have compiled a set of 9 clusters with kSZ measurements from [37l[38l[39] (see table 
[Hand Fig. [3]). Currently both systematic and random errors are very large, but even 
though observations are scarce they already give very interesting bounds on the size of 
a possible void. 

We have used the observations to compare with the constrained GBH model, 
making a grid of models with different parameters (see table [3]) , and calculated the 
likelihood of each model. From the literature we only have the (asymmetric) 68% 
confidence limits for the observations, and not the full probability distributions. To 
respect the asymmetric error bars, while keeping the probability distribution close to 
Gaussian, we use the e-skew- normal (eSN) distribution [10], which has the probability 
density function (PDF) 



f{v;Vi,a ,a+) 



'27r [a 



exp 



a 



V - Vi 



2a 



±2 



(22) 
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where a^ = a~ when v < Vi, and a^ = a^ when v > Vi, to model the individual data 
points. The eSN model is a reasonable model for the data, if we do not have knowledge 
about the detailed likelihood for each data point, since the probability distribution is 
continuous up to the first derivative, the maximum likelihood is at Vi, and 68% of the 
PDF is in the interval [vi — a~,Vi + a^], as required by the data. The mean value is not 
at Vi reflecting the fact that the distribution is skewed. 

When comparing our models to the observations we not only have to take into 
account measurement errors as given in the literature but also the fact that, if the 
model is correct, the apparent peculiar velocities of the clusters are the sums of the 
apparent dipole velocities from the void, and the intrinsic peculiar velocities. In the 
following we assume that the clusters, due to their high redshift, have uncorrelated 
intrinsic velocities, though for future observations of small fields with many clusters this 
may not be true (see e.g. [Illll2lll3]). The peculiar velocities are assumed to be normal 
distributed, and are added in squares to the measurement errors. 

For the full sample (see table [2] for sub samples) the average cluster velocity, 
excluding any systematic shift, and correctly weightecljl is v = 320^400 km s~^, while the 
standard deviation among the cluster velocities is ay = IGSOlggg km s~^. This is a fairly 
large velocity scatter; much larger than expected from linear theory. Given that the 
clusters in current kSZ surveys are among the most massive clusters known, we would 
expect radial peculiar velocities on the order of ~ 400 km s~^ from linear theory |43j . 
which is a factor of 4 lower than what is found from the internal scatter in our data set. 

To estimate the typical peculiar velocity scatter cTpv we will use two different values: 
The average velocity as expected in a ACDM model, fixed at 400 km s~^, and the scatter 
of the current data set, fixed at 1600 km s~^. We could include evolution in the velocities 
— according to linear theory they grow with the growth factor — but since current data 
have very large errors, it would not make much of a difference. 

Summarising, the total log-likelihood for a given model is 

-\ 2 

Vi - av{p, Zi) + tJsys 



21n£ = ^ 



[^.P + <]^/^ 



(23) 



where Vi is the observed value of the velocity, f (p, Zi) is the value according to the model, 
with parameters p, and a^ is a^ ii v > Vi and else a~ . According to [SB] there can be 
significant systematic errors in the data, and we allow for this by adding a systematic 
shift in the values: Wsys- Below we use the systematic shifts t^sys = 0, and t^sys = 750 km 
s~^. The size of the dipole is slightly overestimated by our model. First of all because 
the effect of the void on the CMB sky is not a pure dipole, and second because at large 
distances the projected size of the void on the sky of the cluster becomes less than 27r, 

I We use the eSN distribution for each data point, and weight the data according to Wi = l/cr^, where 
the average standard deviation is the harmonic mean of the error bars: 2/ a = 1/a^ + l/a^. Using these 
weights the PDF of the average velocity is found as fv{v) = ^iUJifi{v)/J2-Wi, while the standard 
deviation is found by Montc-Carlo sampling of the data distributions: ay =< ^^ Wi(vi — u)^/ X^i ^i >j 
with Vi drawn on random from fi{v). 
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Name 


Redshift 


Pec vel 


Syst. Error 


Galactic / 


Galactic h 






[km 


s-i] 




o' 


A1689 37] 


0.18 


+ 1701^^^ 


750 


313.39 


61.10 


A2163 [37] 


0.20 


+490li^° 


750 


6.75 


30.52 


A2261 M] 


0.22 


-l^l^t^f 


750 


55.61 


31.86 


A2390 m 


0.23 


+1900+^i^ 


750 


73.93 


-27.83 


A1835 M] 


0.25 


1 7P-+1675 


750 


340.38 


60.59 


Zw3146 38] 


0.29 


-4nn+37oo 

^UU_i925 


750 


239.39 


47.96 


RXJ1347-1145 39] 


0.45 


1420+}^^° 


? 


324.04 


48.81 


CI 0016+16 3S 


0.55 


-4inn+265o 


750 


201.5 


-27.32 


MS 0451 [38] 


0.55 


4qn+i370 
^yu_ggo 


750 


112.55 


-45.54 



Table 1. Clusters with observed velocities from the literature. The given errors are 
at the l-cr or 68% confidence level. The systematic error is estimated in [3H], and 
is mostly due to confusion with primary CMB anisotropics, and contributions from 
sub-mni point sources. These error sources are not instrument dependent, and we add 
the same error to the observations in [37| and [SHj too. 



Sample 



All clusters 
Cluster in [an [38] 
Clusters in 




^clusters 



9 

8 
6 



Table 2. Statistical properties of sub-samples, and the full sample of clusters. The 
quoted values are in good agreement with what is found in |38| . 



and then the dipole part is even less. To model it we include an empirical rescaling 
factor < a < 1, and to estimate the effect error we use a = 1 and a = 0.8. 

The two choices of each of three free parameters leave us with 8 different likelihoods 
to explore. The two dimensional contours are shown in Fig. HI As can be seen in the 
lower right panel, a void with a central underdensity of Qin = 0.23 (just inside the 3-0" 
limit allowed by other observations, see Fig. [HJand Ref. [S]), cannot be bigger than 1.5 
Gpc at the 3-0" limit. This is the best-case, allowing for generous errors in our modelling, 
an unrealistic large scatter in the peculiar velocities, and a favourable systematic shift 
due to an incomplete analysis of observations. Incidentally, the 3-cr limit for the size 
of the void, allowed by other data, is 1.45 Gpc [19]. For the constrained LTB model, 
interpreting the current kSZ cluster data in the most favourable way, we therefore still 
have at least a 3-0" inconsistency between observations of the geometry of the universe 
(Supernovae, BAO, and 1*^* peak in the CMB), and observations of the kSZ in clusters 
of galaxies. 

On the other hand, if we use the strictest interpretation of the model (upper left 



Looking the void in the eyes - the kSZ effect in LTB models 11 



100 km s"^ Mpc"^ 



0.65 0.44-0.57 0.43 



iiin 



0.13-0.93 



To Ar 

Gpc ro 



0.1-2.5 0.1-0.9 



Table 3. Priors used when scanning the parameters of the GBH models. Hq is a pre 
factor for iifo('') a-nd -ffin and -ffout sue derived from the priors on Qi^ and Hq. 



panel in Fig. H]), then even a void with a shght central underdensity of Qi^ = 0.7 is 
limited to 1.5 Gpc at the 3-cr level, or equivalently a void with a central underdensity 
of fiin < 0.5 have to be less than 1.25 Gpc in size. 

For completeness, in Fig. [5]we also show the equivalent one-dimensional likelihoods. 
However, we caution against overinterpreting them: current data gives a non-detection 
of a local void with four clusters moving towards us and five moving away, and the 
two-dimensional likelihood contours in the tq-Q-u^ plane are strongly degenerate. Hence, 
the reduced one-dimensional contours give unreasonably large confidence limits, and 
we refrain from quoting any general bounds on single parameters, but conclude that 
measurements of the kSZ effect in clusters is by far the most constraining way to limit 
the size of a local void described by an LTB metric. 

5. Future Experiments 

In this section we explore the limits on the LTB models placed by future experiments 
like the South Pole Telescope [S], the Atacama Cosmology Telescope |1S], the APEX 
telescope [l6], or the Planck satellite [17]. Below we use parameters for the ACT, but 
our results apply equally well to any of the four surveys. Our basic hypothesis in this 
section is that the Universe is homogeneous, and we show that giga parsec sized voids 
with a significant underdensity are strongly ruled out, even with a few observed galaxy 
clusters. 

The first generation of large scale kSZ survey is characterised by a limited number 
of frequency bands (ACT, APEX and SPT), or by low angular resolution (Planck), 
and even though the resulting cluster catalogues will be a gold mine for studying the 
thermal SZ signal, it is much harder to exploit the full potential of the kinematic SZ 
for a number of reasons: the intensity dependence is similar to that of the primordial 
CMB, point sources can contaminate the fluxes, giving a need for excellent resolution, 
and finally what is really measured is the relative dip (or bump) in the frequency due 
to the integrated (Thomson) scattering of photons. That is, what is measured directly 
is the momentum of the cluster gas with respect to the CMB, not the peculiar velocity 
itself, see Eq. (I19p . If only three (ACT) or four (SPT) frequency bands are available 
then an estimate of the total gas mass, and the internal velocity dispersion, are needed 
in order to extract the peculiar velocity. The gas density can be estimated using X-rays 
together with thermal SZ measurements, or alternatively using weak lensing together 
with extensive modelling, while many precise redshifts of cluster galaxies in each cluster 



Looking the void m the eyes - the kSZ effect m LTB models 



12 




Op,=400,V3,,=0,a=0.8 : 




z.i> 
2.0 




> n 


o 1.5 
(5 


- ^,^ 


C^ - 


1.0 


. A-O 


"^^"^^ - 


0.5 


- 


ap„=1600,V3y3=0,a=0.8: 



0.2 



0.4 0.6 



0.8 



0.2 



0.4 0.6 





Figure 4. Likelihood contours for current observed clusters with 1-ct to 3-(t regions 
given by bright to dark yellow. From left to right, top to bottom we have cTpv = 
(400, 1600) km s-\ Usy^ ^ (0, -750) km s~\ and a = (1, 0.8). 
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Figure 5. One dimensional likelihoods for current observed clusters with 1-cr and 2-cr 
levels indicated by dotted lines. Velocities are in km/s. 



are needed to constrain the internal velocity dispersion [l8] . The bottom line is that even 
though 10^-10^ clusters will be detected in the first generation of large scale SZ surveys 
we can only constrain the peculiar velocities of those clusters that are well measured 
by other means, severely limiting the clusters that can be explored, to those at lower 
redshifts {z < 0.4). Hence the number of clusters with kSZ measurements will be close 
to the number of clusters which are well observed in X-rays, i.e. 10^-10'^ clusters. 
To construct a mock survey we model the cluster PDF as 



n{z) oc 



exp 



5/3' 



exp 



(- 



(24) 



in accordance with the cluster density given in [S] if we set 2max = 0.6, see Fig. El 
The last term is an exponential redshift cut-off, due to the limited reach of X-ray 
observations, and we set zx = 0.4. Notice that the detailed form of the cluster PDF 
is not so important because X-ray observations mostly sample redshifts lower than the 
maximum where the volume density is nearly constant, and the redshift density then 
goes like z^. 

We use a similar likelihood function for our mock survey as the one in Eq. 



2\nC = J2 



Vi - av{p, Zi) + V, 



sys 



\(J- 



a2 1V2 

pvJ 



(25) 



where the only difference is that we assume Gaussian errors for our observations. Due 
to uncertainties in the internal velocity scatter in the cluster gas, confusion with the 
primordial CMB and confusion with point sources in any future experiment it will be 
very difficult to reduce the uncertainty in the peculiar velocity Uj to less than ~ 100 km 
s~^ [IH]. In the first generation experiments it will be hard to get an uncertainty ctj < 500 
km s~^, due to the limited number of frequency bands, uncertainties in measuring the 
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Figure 6. Cluster redshift number density used to select the redshift distribution 
for forecasting future surveys. The dashed line is the number density folded with the 
X-ray cut-off. 



gas mass by other means, and incomplete knowledge of point sources in the clusters 
Below we will consider a best-case scenario of ai = 400 km s~^, and a worst case of 
800 km s~^. For simplicity we will fix the rescaling of the modelled velocities due 



a. 



to uncertainties in the theoretical model to a = 0.8, and we will consider systematic 
errors of Ugys = and 400 km s~^, while the scatter in the peculiar velocities is fixed 
as above to CTpv = 400 and 1600 km s~^. These are reasonable choices that covers the 
probable spread in both the observable and the theoretically motivated parameters in 

Eq. (I25D. 

We use Monte Carlo modelling to sample the PDF of the cluster redshift distribution 
Eq. (121]), and we assume that the clusters have normal distributed peculiar velocities 
Vi ~ A/'(0, (Tpv). In Figs.[7]and[8]are shown the 2D likelihoods in the cases were we have 10 
and 100 clusters with good kSZ measurements. As can be seen already with 10 clusters, 
strong bounds can be put on the LTB model, while with 100 clusters, assuming the FRW 
model is correct, giga parsec sized voids are essentially ruled out. Complement arily we 
have shown that, if the measurements have systematic errors at the 400 km s~^ level, a 
spurious positive detection of a very shallow underdensity (fijn < 1) could be inferred. 
Conversely, care should be taken when constraining LTB models with kSZ observations 
where systematic errors have been included self-consistently in the analysis, under the 
assumption of a homogeneous universe, by fixing the average peculiar velocity to zero, 
as advocated in Ref. 



6. Discussion and conclusions 



As an alternative explanation of the dimming of distant supernovae it has recently been 
advocated that we live in a special place in the Universe near the centre of a large 
void. The universe is no longer homogeneous and isotropic and the apparent late time 
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Figure 7. Likelihood contours for a mock survey witli 10 well observed clusters with 
1-cr to 3-cr regions given by bright to dark yellow. Also indicated in the individual 
plots are the parameters used for the simulation. Velocities are given in km s~^. 
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Figure 9. Likelihoods for the constrained model for complementary data. The 
combined data sets of Supernovae Type la (SNIa), BAO and WMAP arc shown in 
yellow with 1-. 2-, and S-cr contours, while the individual SNIa, BAO, and CMB data 
sets are shown in blue, purple, and green respectively with 1- and 2-cr contours. The 
figure is taken from [19]. 



acceleration is actually a consequence of spatial gradients in the metric. If we did not 
live close to the centre of the void, we would have observed a CMB dipole much larger 
than that allowed by observations. Hence, until now it has been argued, for the model 
to be consistent with observations, that by coincidence we happen to live very close to 
the centre of the void. However, even if we are at the centre of the void, we can observe 
distant galaxy clusters, which are off-centre. In their frame of reference there should 
be a large CMB dipole, which manifests itself observationally for us as a kinematic 
Sunyaev-Zeldovich effect. 

In this paper we have studied the induced dipole for off-centred clusters due to the 
different trajectories of photons from the last scattering surface, and computed the size 
of the corresponding apparent velocity of those clusters with respect to the rest frame 
of the CMB LSS, depending on the parameters of the void model. We then analysed 
the present observations of the kSZ effect in a handful of clusters and gave very strong 
constraints on the size of the void in LTB models. In fact, for our specific constrained- 
GBH model the bounds are in conflict with other constraints from supernovae, baryon 
acoustic oscillations and CMB, at least at the S-a level, see Fig. |9l and therefore we 
conclude that the constrained-GBH void model is practically ruled out if the current 
interpretation of kSZ observations are correct. 

At present, kSZ data leave small voids (rg < 800 Mpc) unconstrained, 
independently of the inner density contrast (fiin) used, simply because the radius is so 
small that the void does not impact the cluster with the lowest redshift in the sample, 
see Table 1. In order to put limits on small voids with large density contrasts, or sudden 
transitions, we would need to use local large scale structure data that can be checked 
for homogeneity, or alternatively low redshift supernovae. 

Current kSZ observations are limited in numbers, and are still not precise enough 
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to make a single positive detection of the peculiar velocity of a cluster at the 2-cr level, 
see Fig. El Hence, one cannot rule out that they are plagued by large systematic 
and/or random errors. This is about to change in the coming years, and we have made 
predictions of how systematic near-future observations from kSZ surveys like ACT or 
SPT could strongly rule out all LTB models with giga parsec sized voids. In the case 
that the LTB models were confirmed, the average apparent velocity profile as a function 
of distance would give a direct handle on the density and expansion rate of the void. 
This unique relationship makes by far kSZ observations the most powerful data for 
constraining LTB models. 

If the sensitivity of experiments like ACT, SPT and Planck are as planned, and the 
present systematic errors are under control, it is expected that large voids will definitely 
be ruled out at many sigma. If these experiments yield 10 (100) well observed clusters 
we could reasonably expect to rule out voids of 800 (500) Mpc radius at the 3-0" level. 
We hope to come back to this severe challenge to LTB models and redo the analysis 
with the new data in the near future. 
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